Realizable spin models and entanglement dynamics in superconducting flux qubit systems 
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Realizable spin models are investigated in a two superconducting flux qubit system. It is shown that a specific 
adjustment of system parameters in the two flux qubit system makes it possible to realize an artificial two-spin 
system that cannot be found naturally. For the artificial two-spin systems, time evolution of a prepared quantum 
state is discussed to quantify quantum entanglement dynamics. The concurrence and fidelity as a function of 
time are shown to reveal a characteristic entanglement dynamics of the artificial spin systems. It is found that the 
unentangled input state can evolute to be a maximally entangled output state periodically due to the exchange 
interactions induced by two-qubit flipping tunneling processes while single-qubit flipping tunneling processes 
plays a role of magnetic fields for the artificial spins. 
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Introduction. Superconducting qubit systems as one of 
promising candidates have been paid much attentions for 
quantum information processing and computing. The tunable 
superconducting devices have provided a variety of possibili- 
ties to realize quantum spin models that are not findable nat- 
urally. Recent experiments have shown that different types of 
exchange interactions are observable. Particularly, there have 
been demonstrated an Ising type interaction in two charge 
qubits HJ] and two flux qubits J 2, 13] and an XY type inter- 
action in two phase qubits 0, □ ] as well as superconducting 
single qubits iaullSl- Moreover, such realizations of artifi- 
cial spin systems makepossible to observe entangled states 
of two qubits |Q1 0, B Si ■ Indeed, for the time evolution of 
states in the experiments of charge [1] and phase qubits |4|], 
a partial entanglement has been observed. An experiment of 
a capacitively coupled two phase qubits |5] shows that higher 
fidelity for the entanglement exhibits in an excited level. The 
higher fidelity is caused by two-qubit tunneling processes [9] 
between two qubit states, i.e., flipping both qubits. Such a 
two-qubit tunneling processes contributes exchange interac- 
tions between the two artificial spins. 

In this paper, we will theoretically investigate a possible 
realization of quantum spin models in superconducting flux 
qubit systems by varying a system parameter. Especially, we 
use a phase coupling by introducing a connecting wire be- 
tween the two qubit loops (see Fig. Q} | flol] because the phase 
coupling gives more controllable parameters than the induc- 
tive coupling with respect for the manipulation of qubit states 
Hill . It is shown that the Josephson junction in the connecting 
superconducting wires plays a role of controller in determin- 
ing exchange interactions between the two qubits. In general, 
it is found that the two flux qubit system can map into an XYZ 
quantum spin model in the presence of magnetic fields. We 
show that specific values of system parameters generate vari- 
ous types of quantum spin models. Further, to address about 
time evolution of an input state for the two flux qubit sys- 
tem corresponding to quantum spin models, we introduce the 
concurrence and fidelity as a function of time as a measure 



of entanglement and evolution of the state. It turns out that an 
unentangled (entangled) input state evolves to be an entangled 
(unentangled) state periodically with a characteristic period of 
time. 

Time evolution of quantum states. A system described by 
two quantum states can be a qubit. The two states can be rep- 
resented in terms of pseudo-spin language, i.e., two orthogo- 
nal states |t) and |X). Then, any normalized pure state of two 
qubit systems can be written as a linear combination in the 
basis {|TT), IU>, lit), IU»: 



\^) = a\n)+b\n)+m)+d\u). 



(i) 



For two qubit systems, a given Hamiltonian H generates the 
time evolution of the state through the Schrodinger equa- 
tion ihd t \ifr(t)) = H \4i{t)). If the Hamiltonian H is indepen- 
dent of time, the time-dependent state is given by \tf/(t)) - 
exp [— ^ ] l<K0))> where \ifr(0)) = \ifr) is a given state at the ini- 
tial time. By virtue of the unitary transformation U making 
the Hamiltonian diagonal, at time f, the state is given by 



\m)=G(t)mo)), 



(2) 



where the propagator is G(t) = U exp [~ ^p^ ] for me 
time evolution of the state. 

To quantify entanglement for the time evolution of the state, 
we introduce the concurrence as the overlap between the state 
and the spin flipped state at a given time t : 



c(i<Mf)» = \(m i m)\ ., 



(3) 



where the spin flipped state is given by \\p{t)) = cr y ® cr y \iff*(f)) 
lfl2ll with the pauli matrix o~ y . The concurrence ranges from 
zero (unentangled state) to one ( a maximally entangled state). 
To help understanding the entanglement dynamics, one can 
define the overlap between the states at the initial time (input 
state) and at a given time t (output state) as the fidelity: 



Fit) = \(m i <ko)>|. 



(4) 
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FIG. 1 : (color online) A two flux qubit system. The system is com- 
posed of two (left and right) qubit loops. In order to couple the two 
flux qubits, we use two connecting superconducting wires where the 
Josephson junction E'j plays the important role for controlling the in- 
teraction between the two qubits since the two wires give the bound- 
ary condition as a function of phases {<p^,ipi, tp'} from the fluxoid 
quantization along the closed path through the two connecting wires. 
By varying the amplitude of E'j, the two flux qubit sytem can be map 
into a quantum two-spin model. The state of each qubit loop is in a 
superposed state of which |J,) and |f ) represent the diamagnetic and 
paramagnetic current states, respectively. This schematic of the sys- 
tem show the state |TT) that is one of possible four states. Here, O 
and ® denote the direction of the magnetic fields, /ip) = ©lp)/^. m 
the qubit loops. Eji, Ej, E'j are the Josephson coupling energies of 
the Josephson junctions in the qubit loops and the superconducting 
connecting wire, and ip's are phase differences across the Josephson 
junctions. 



If F(T) = 1, the output quantum state is the same with the 
initial input state at t = T, i.e., unentangled (entangled) initial 
state returns to unentangled (entangled) state. Then, for time 
evolution of quantum states, entanglement dynamics can be 
understood from the concurrence and fidelity. 

Model - We consider a two superconduting flux quit sys- 
tem where the two flux qubits interact each other by control- 
ling the Josephson junction energy E'j in the superconducting 
connecting wire in Fig. [TJ We assume that the inductances 
associated with the geometry of the system is so small that the 
inductive energy is negligible. The Hamiltonian describing 
the model is given by the sum of the charging and Josephson 
energies: 

H({<pi, Cp\ <fii, ip'}) = H c ({<Pi, ?'}) + Hj({<p u <p'}), (5) 

where the phases across the Josephson junctions are ipi and 
their time derivatives are <pi. The charging energy of Joseph- 
son junctions in the two qubit loops and the connecting wire 
is given by 



x ' V i=l ae{a,b,c) 



(6) 



where C a (C) are the capacitance of the Josephson junctions 
in the qubit (connecting) loops. (Do = h/2e is the unit flux 
quantum. The Josephson energy of the junctions is given by 



= Z Z 2£>n 2 ^ + 2£;sin 2 ^, 



(7) 



/'=! a€{a,b,c] 



where E/s are the Josephson energy of junctions in the qubit 
and connecting loops. For flux qubits, the charging energy 



is much smaller than the Josephson energy. The number of 
Cooper pairs n and the phase tp are non-commuting variables, 
i.e., [if, n] = i, such that the canonical momentum can be 
introduced as P„ = n% = -ifid^, where n - q/2e with the 
charge from the Josephson relation q = C(<J>o/2^)(^. At low 
energies, then, the charging energy Hc({<fiY) play a role of ki- 
netic energy for a phase particle, while the Josephson energy 
Hj({(fii}) plays a role of confinement potential for the particle. 

From the fluxoid quantization along the closed paths in the 
two qubit loops and the connecting loop path including the 
two superconducting wires, the constraint conditions for the 
phases are given by 



ae\a,b,c] 

2nn' + ((p$ - ipl) - ip' = 0, 



(8) 
(9) 



where n, and n' are an integer and f, = <t)/<l>o are the applied 
flux in the qubit loop i. The constraint of tp" and ip a 2 in Eq. (0 
induce the coupling between the two flux qubits, which is call 
the phase coupling. 

In the low energy limit, generally the Hamiltonian of su- 
perconducting flux qubit systems can be written in terms of 
the circulating current states in each qubit loop 11311 . For two 
flux qubit systems, following Ref. [9, , one can write the 
two-qubit matrix Hamiltonian in terms of qubit energy levels, 
single-qubit tunnelings, and two-qubit tunnelings: 



f E 



H = 



TT _ 'i 



f 2 E n 



-d -t\ -t\ £ 



-h 

II ) 



(10) 



where £"s are the energies for the two qubit states. 



qubit states correspond to the local minima \w° 

and ipf c 
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Josephson energy Hj(tpi) with E' 

Here, m, =T and { for the qubit i. The two-qubit level en- 



ergies are given by E„ 



,})• 



where the characteristic oscillating frequencies are cof. mi m , = 

W t ^ H ^M^ with an effective mass M < = ©Nff 
and effective capacitance in the harmonic oscillator ap- 
proximation 11411 . t\ and ti are the single- and two-qubit tun- 
nelings between the two states of two qubits. Single-qubit 
tunneling describes single-qubit flipping for the macroscopic 
quantum tunneling between the two states of the two qubit 
states. For example, |fT) \\Vl- The two-qubit tunnel- 

ing amplitudes, (i) and (ii) f*, describe the tunneling pro- 



cesses, (i) |TT) <==> in the parallel pseudo-spin states and 
(ii) |T4) <==> IIT) m the anti-parallel pseudo-spin states. The 
tunneling amplitudes are calculated by the numerical methods 
such as WKB approximation, instanton method, and Fourier 
grid Hamiltonian method 1 15]. 

In fact, the tunneling amplitudes and the low energy qubit 
energies are determined by the systems parameters of the su- 
perconducting flux qubit system. Once the parameters are ad- 
justed, generally, an artificial spin Hamiltonian is given in a 
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form from Eq. ( [Tot : 



2 (11) 



ae{.ic,y,z| 



where = -fj 

(•Err 



B 



En+E 



-Ex-E u )/4, 
, r -%)/4, 7, = -(^ + <*)/2, J y = 
% + E ll) / 4 - The sin 8 le 



{f 2 - f*)/2, and 7, = (E n - E- 
qubit tunnelings play the role of a transverse magnetic field 
while the energy difference of two-qubit levels correspond to 
the applied magnetic field parallel to the z-direction of spins. 
Note that the x- and y-components of the exchange interac- 
tion are determined by the two-qubit tunnelings and the z- 
component of the interaction is the energy difference between 
the parallel spin state and the anti-parallel spin state. Conse- 
quently, Eq. QT]shows that an XYZ quantum spin model with 
magnetic fields can be realizable in any two flux qubit system 
111 311 . We will discuss a specific realization of a quantum spin 
model with adjusted system parameters and entanglement dy- 
namics for an input state. 

Realizable artificial spin systems. Case I. For E'j — O.OEj 
and Eji — 0.7 Ej, a two-spin Hamiltonian can be constructed 
by the relations of E^ = = = E^ and ^ = ^ = ti. 
The numerical values of the macroscopic quantum tunnelings 

0.00024£ / . From 



Mb) 



are obtained as t\ = 0.0075£y and t-~ ' = j. 
Eq. QT| the two flux qubit system is described by the corre- 
sponding spin Hamiltonian: 



(12) 



where B — —t\ — and J — —t 2 . Note that the entanglement 
dynamics of this spin system is determined only by the single- 
and two-qubit tunneling amplitudes. 

For the spin system, the concurrence is given by 



C(t) = [C + Cicos47f] 



1/2 



(13) 



where Co = [(a + d) 2 -(b + c) 2 ] 2 /4 + [(a-d) 2 -(b-c) 2 ] 2 /4and 
d = -[(a + d) 2 -(b + c) 2 ][(a-d) 2 -(b-c) 2 ]/2. It should be 
noticed that the concurrence does not depend on the magnetic 
field B = —t\, i.e., the single-qubit tunneling. The concurrence 
is an oscillating function with respect of the exchange interac- 
tion / = -?2 with the characteristic period of time T = n/2J. 
This shows that any entangled state cannot be generated by 
applying the magnetic field in this quantum spin system. At 
t = (m + 1)tt/2J with an integer m, the concurrence reaches its 
maximum value C(f) = 2\ad — bc\, i.e., a maximally entangled 
state. At t — ^jcos" 1 the entanglement disappears, i.e., 
the input state involves to be unentangled. 

The fidelity of this quantum spin system is given by 



F + V F'[ cos 2(5 + crJ)t + F 2 cos 4Bt 



1/2 



(14) 



where F = 1 - [(a + d) 2 + (b + c) 2 ] [(a - d) 2 + (b - c) 2 ] /2 - [(a + 
df-(b + c) 2 ] 2 l%,F\ = (a+b + c + d) 2 [(a- d) 2 + (b - c) 2 ]/4, 



F~ = (a - b - c + d) 2 [(a - d) 2 + (b - c) 2 ]/4, and F 2 = [(a + 
d) 2 - (b + c) 2 ] 2 /8. For zero magnetic field, the quantum state 
evolves in time due to the exchange interaction. The fidelity 
has twice longer period of time than the concurrence. 

Let us study entanglement when the initial input states are 
not in an entangled state. We choose the case of c = d = 0, 
i.e., PF(0)> = a |f T) + b |f |>, in which the initial state is a 
product state. It is clearly shown that the exchange interaction 
makes the artificial spins entangled in the concurrence C(f) = 
|(2a 2 - 1) sin2/f|. The period of the concurrence is T — n/J. 

For a — d and b — c, the initial state can be written in 
the Bell basis with = T=(|f l> + lif » and ^ = 

i(|ff> ± |||». The state is given by |<A(0> = V2(a |0 + > + 
b | V F + )). Initially the input state is an entangled state quantified 
by the concurrence C(0) = 2|a 2 - b 2 \. For the time evolution 
of the state, the concurrence is not a function of time, i.e., 
C(?) = 2|a 2 - b 2 \. However, the fidelity is a function of time. 
It does not depend on the exchange interaction, i.e., F(t) = 
VA + A x cos4Bf, where A = 1 - 2(a 2 - b 2 ) 2 and A\ = 2(a 2 - 
b 2 ) 2 . This can be understood as follows. The initial state can 
be rewritten in the eigen basis: \if/(t)) = (a + b) \i//q) + (a - 
b) |i^3), where \if/o) and 1^3) are the ground state and the third 
excited state, respectively. At a given time t, then, \if/(t)) = 
(a + b)e~' E "' \iffo) + {a- b)e~' Ei ' where the energies for the 
ground state and he third excited state are £0 = J - 2B and 
£3 = J + 2B. As a result, for the constant concurrence, the 
fidelity is oscillating in time. 

Case II. For E'j = Q.6Ej and Ej\ = 0.7 Ej, one find the 
relations: E^ = En and E^ - E^, and t\ = = 0. The 
corresponding two spin system can be written as 



H = J(S*Si 



S j^) ~\~ J r- S 1 5"^. 



(15) 



where J = —t%/2, and J z = (E^ - En)/2. The numerical val- 
ues of the two-quit tunneling amplitude and the energy differ- 
ence between the two states are given by = 0. 00024/? j and 
J z = -0.425045 E j. The Hamiltonian describe an anisotropic 
spin exchange interaction belonging in the class of the XYZ 
spin model. Interestingly, the x- and z-components of the in- 
teraction are anti-ferromagnetic because / < and J z < 
while the y-component is ferromagnetic. To our knowledge, it 
is unlikely to find a class of spin Hamiltonian naturally. 
For the spin Hamiltonian, the concurrence is given by 



at) 



Co + Y cos 4(7 + a-J z )t + C 2 cos SJt 



1/2 



(16) 



where Co = 
Cj = 2(a 
given by 



[(a + d) 4 + (a-d) 4 ]/4 + 4b 2 c 2 ,Ct = -2{a + dfbc, 
- dfbc, and C 2 = -(a 2 - d 2 ) 2 /2. The fidelity is 



F(t) 



F + ^ Fi cos 2(7 + crj z )t + F 3 cos4/f 



1/2 



(17) 



where F = 1 - 2(a 2 + d 2 )(b 2 + c 2 ) - (a 2 - d 2 ) 2 /2, F+ = (a + 
d) 2 (b 2 +c 2 ), F~ = (a-d) 2 (b 2 +c 2 ), andF 2 = (a 2 -d 2 ) 2 /2. This 
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FIG. 2: (color online) Time evolutions of the concurrence C(t,9) 
and the fidelity F(t,6) for the input state ^(O)) = cos2;r0|TT) + 
sin 27T0 14.1) at the co-resonance point f\ = fi = 0.5 in the two su- 
perconducting flux quit system, (a) Case I. For E'j - Q.QEj and 
Ej[ = O.lEj, the two flux qubit system corresponds to the spin 
Hamiltonian H = J S\ S* + B (S\ + Sj) with the exchange in- 
teraction J = -f, and the magnetic field B = —t\. The single- 
and two-qubit tunneling amplitudes are given by t t = 0.0075£/ and 
t2 = 0.00024£/. The characteristic period of time is T\ = n/2J for 
the both concurrence and fidelity, (b) Case II. For E'j = 0.6Ej and 
Ej[ = Q.lEj, the two flux qubit system maps into the spin Hamilto- 
nian H = J (5 J SS- S\ S y 2 ) + J z S\ S 2 with the exchange interaction 
J = -t%/2 and J, = (E^ — £ T jJ/2. The two-quit tunneling amplitude 
is f 2 ' = 0.00024.E/ and the energy difference between the two states 
is J- = -0.425045£/ . The period of time is T 2 = n/2J. The fidelity 
has twice period of the concurrence, (c) Case III. For E'j = 0.05Ej 
and Ej t = Q.lEj, The two flux qubit system is described by the spin 
Hamiltonian H = J{ S\ 5* - S\ S y 2 ) + J z S\ S\ + B (S\ + S £) with the 
magnetic field B = —t\ and the exchange interactions J = -t^/2, and 
J z = (£ff - Efi)/2. The single- and two-qubit tunneling amplitudes 
are ti = 0.0024£/ and f 2 = 0.00024£/ and the energy difference be- 
comes J- = -0.05Ej. The characteristic period of the time evolutions 
is T3 = 0.687r/27. The fidelity has twice period of the concurrence. 



shows that the concurrence and fidelity have a similar dynamic 
property. However, the fidelity has twice longer period than 
the concurrence. 

Compared to the Case I, for the initial product state |*P(0)} = 
a ITT) + b\1i}, the concurrence becomes C(f) = a 2 |sin47f| 
with the period T = n/2J. The input state \iff(t)) = 
V2(a |<1> + ) + b\ x i l+ }) evolutes and its concurrence is oscillat- 
ing with C(f) = 2 V« 4 + b 4 - 2a 2 b 2 cos 4(7 - J z )t . Also, the 
fidelity is given by F(f) = V^o + A i cos 4Jt , where Aq = 



1 - 8a 2 b 2 and A\ = %a 2 b 2 . Thus, the Case I and II show a 
different characteristic entanglement dynamics depending on 
the realizable two artificial spin models in the flux qubit sys- 
tems. 

Case HI. For E'j = 0.05£/ and E n = O.lEj, The two- 
qubit energies have the relations: Eff = En, and E^ = Ef±. 
The two-qubit tunneling becomes A = 0. We find another 
realization of a two spin system: 

H = J (S? S x 2 - S\ S\) + J z S\ S\ + B (S I + S 2 C ), (18) 

where B = -t\, J = -f°/2, and J, = (E n - E n )/2. For 
the system parameters of the superconducting flux qubits, the 
single- and two-qubit tunneling amplitudes are t\ = 0.0024£' j 
and t% = 0.00024£'y and the energy difference is J z = 
-Q.Q5E j. The expressions of the concurrence and fidelity are 
too lengthy to display. 

In Fig. [2] we plot the concurrences and fidelities as a func- 
tion of time t and the initial state parameter 6 to give the com- 
parison of entanglement dynamics between the three different 
spin models for the same initial state | V P(0)> = cos27r6*|TT) + 
sin 2n6 \H). Explicitly, the different values of system param- 
eters controlling the two flux qubits are given in the captions 
of the figures. For the time evolution of the initial state, it is 
shown that unentangled (entangled) state can become entan- 
gled (unentangled) state even though the specification of the 
superconducting devices are different each other. 

Summary. A two superconducting flux qubit system has 
been considered to investigate a possible realization of quan- 
tum spin models. Three different artificial spin models were 
demonstrated by varying controllable system parameters. The 
realizable spin models in the flux qubit system are not likely 
to find naturally. We discussed the entanglement dynamics 
of the artificial spin models in the specific parameter values 
of the two superconducting flux qubit system. It was found 
that the input unentangled (entangled) state can become max- 
imally entangled (unentangled) state irrespective of the spec- 
ifications of the superconducting devices. Such a maximally 
entangled state should be observable experimentally. 

We thank Huan-Qiang Zhou and John Paul Barjaktarevic 
for helpful discussions. 
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